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Abstract. - The method of detection of the unstable periodic spatio-temporal states of spatial 
extended chaotic systems has been proposed. The application of this method is illustrated by 
the consideration of two different systems: (i) the fluid model of Pierce diode being one of the 
fundamental system of the physics of plasmas and microwave electronics and (ii) the complex 
one- dimensional Ginzburg-Landau equation demonstrating different regimes of spatio-temporal 
chaos. 




It is well known that the unstable periodic orbits 



^—j- (TJPOs) embedded into chaotic attractors play an impor- 
^f*} tant role in the dynamics of the systems with a small num- 
<\j~ber of the degree of freedom [1-3]. The chaotic regime of 
OQthe system may be characterized by means of the set of 
<^> UPOs [4]. A universal and powerful tool for exploration of 
C**- chaotic dynamics [5] , unstable periodic orbits proved to be 
especially efficient in the context of chaotic synchroniza- 
£>tion [6-9]. The different types of chaotic synchronization 
J^j (such as phase synchronization [8], lag [10] and complete 
^Jj synchronization [11]) may be explained in terms of unsta- 
C^ble periodic orbits [12]. Eventually, UPOs play the key 
role for the chaos controlling problem [13] since the unsta- 
ble periodic orbits may be stabilized by means of the week 
influence on the system dynamics, e.g., with the help of 
small variation of the control parameter [14] or with the 
feedback of different types [15-17]. 

In the spatial extended systems the unstable periodic 
spatio-temporal states (UPSTSs) exist [18] which are sim- 
ilar to the unstable periodic orbits in the chaotic systems 
with a small number of the degree of freedom. In particu- 
lar, the chaotic dynamics of spatial extended systems may 
be controlled by stabilizing such unstable periodic spatio- 
temporal states [19]. Therefore, one of the important 
problem connected with the study of the spatial extended 
chaotic system is finding these unstable periodic states. It 
is appropriate to suggest that the methods aimed at the 
search of UPOs of discrete maps (and the flow dynamical 
systems with small dimension of phase space, too) may 
be adapted to the spatial extended systems. The method 



proposed by D.P. Lathrop and E.J. Kostelich [2], as an ex- 
ample, had been used to pick out UPSTSs for the fluid 
model of Pierce diode [20]. This method is based on the 
obtaining the histograms describing the frequency of sys- 
tem returning to the vicinity of UPOs (in the systems with 
a small number of the degree of freedom) or UPSTSs (in 
the spatial extended systems), respectively. Nevertheless, 
this method applied to spatial extended systems is rather 
imprecise and time-consuming. Let us also note the work 
of S.M. Zoldi and H.S. Greenside [21], where the numeri- 
cal analysis of UPOs for a high- fractal-dimension chaotic 
solution of the partial differential equation is carried out 
with the help of the innovative damped-Newton method. 

In this Letter we describe the modification of the 
method of P. Schmelcher and F. Diakonos (SD-method) 
[22,23] allowing precise detection of UPSTSs in the spatial 
extended chaotic systems. As the sample analyzed spa- 
tially extended chaotic systems we consider here the fluid 
model of Pierce diode and the complex one-dimensional 
Ginzburg-Landau equation. 

As the primary system under study we have used the 
fluid model of Pierce diode [24-29] being one of the 
simplest beam-plasma systems demonstrating complex 
chaotic dynamics. It consists of two plane parallel infi- 
nite grids pierced by the monoenergetic (at the entrance) 
electron beam. The grids are grounded and the distance 
between them is L. The entrance space charge density 
po and velocity vo are maintained constant. The space 
between the grids is evenly filled by the neutralizing ions 
with density \pi/po\ = 1- The dynamics of this system is 
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defined by the only parameter, the so-called Pierce param- 
eter a = ujpL/vo, where uo p is the plasma frequency of the 
electron beam. With a > tt in the system, the so-called 
Pierce instability [27, 30] develops, which leads to the ap- 
pearance of the virtual cathode. At the same time, with 
a ~ 37r, the instability is limited by non-linearity and the 
regime of complete passing of the electron beam through 
the diode space can be observed. In this case the system 
can be described by the partial differential equations: 



dv dv dip dp dp dv ^ 
dt V dx dx' dt V dx ^ dx 



5V 
dx 2 



a 2 (p-l), 



(1) 
(2) 



with the boundary conditions: 



v(0, t) = 1, p(0, t) = 1, y>(0, t) = y>(l, t) = 0. (3) 

In equations (1) the non-dimensional variables (space 
charge potential ip, density p, velocity v, space coordinate 
x and time t) are used. They are related to the corre- 
sponding dimensional variables as follows: 



iff = (v$/r,)<p, E> = (v*/L V )E, 
pf = pop, v' = v v, x' = Lx, t' = (L/v )t, 



(4) 



where the dotted symbols correspond to the dimensional 
values, r] is the specific electron charge, an d po are the 
non-perturbed velocity and density of the electron beam, 
L is the length of the diode space. Equations (1) are inte- 
grated numerically with the help of the one-step explicit 
two-level scheme with upstream differences and Poisson 
equation (2) is solved by the method of the error vector 
propagation. The time and space integration steps have 
been taken as At = 0.003 and Ax = 0.005, respectively. 

One of the core problems related to the spatial ex- 
tended system consideration is the infinite dimension of 
the "phase space" W°°. As a consequence, the state 
XJ(x,t) of the system of investigation should be considered 
instead of vector x(t) in R n as in the case of the flow sys- 
tems 1 . After the transient finished (i.e., t > t tr ) the set of 
the states U Vt > t tr may be considered as attract- 
ing subspace W s of the infinite-dimensional "phase space" 
W°° of the spatial extended system under study. If the 
dimension of this subspace is finite, the finite-dimensional 
space R m of variables may be used to describe the dynam- 
ics of the spatial extended system. 

In is well-known that SD-method was developed to the 
UPOs detection in the systems with discrete time, al- 
though it may be also applied to the flow systems [23] 
by means of reducing them to maps with the help of 
Poincare secant. In order to apply the SD method to an 
extended system, we assume that its infinite-dimensional 
phase space possesses the low-dimensional attracting in- 
variant subspace W s , and the desired solution lies in this 



1 For the system under consideration (1) this state XJ(x,t) is the 
vector of the functions characterizing the system dynamics, i.e., 
U(x, t) = (v(x, £), p(x, t), <p(x, t)) T . 



subspace. Further, we construct the auxiliary system y(t) 
in which the vector field y is in one-to-one correspondence 
with W s . 

The stationary states XJ°(x,t) = XJ°(x) of the spa- 
tial extended system correspond to the fixed points in 
the phase space of the auxiliary system, while the pe- 
riodic spatio-temporal states of (1), (2) are in one-to- 
one correspondence with the periodic orbits of the finite- 
dimensional system y(t). Therefore, UPSTSs of spatial 
extended system may be found by means of the detection 
of UPOs of the auxiliary finite-dimensional system. 

There are many well-known methods for applying low- 
dimensional variable space to describe the behavior of the 
spatial extended system, among which a typical one is 
the mode expansion method. In particular, in Ref. [28] 
the low-dimensional model has been constructed for Pierce 
diode (l)-(2) by means of the extraction of several prin- 
cipal modes with the help of Galerkin method. In the 
present work we propose the use of the variables taken 
from several points Xi of the extended system space to 
construct the finite dimensional system 



y{t) = (p(x 1 ,t),...,p(x m ,t))' 2 



(5) 



where m is the dimension of the auxiliary system, 
Xi = iL/{m + 1), i = 1, m. In comparison with the other 
known methods, such approach allows us to undergo eas- 
ily from the spatial extended system state U(x, t) to the 
low-dimensional vector y(t) without any additional calcu- 
lations or measurements. 

For the system under study (l)-(2) we have estimated 
the dimension of the auxiliary vector y(t) as m = 3. This 
assumption is based on the previous results of the consid- 
eration of the finite-dimensional model of the Pierce diode 
dynamics obtained with the help of Galerkin method [28]. 

To confirm meeting of the requirements of the one-to- 
one correspondence between state XJ(x,i) of the spatial 
extended system and vector y(t) of the constructed aux- 
iliary system with the small number of degree of freedom 
we have used the neighbour method [31]. We have exam- 
ined that the distance d(y± , y2) = | |yi — y2 1 1 between two 
vectors yi = y(ti) and y2 = y(^) taken in the arbitrary 
moments of time t\ and t<i is close to zero if and only if the 
distance 5(Ui,U2) between two different states U(x,ti) 
and U(x,t2) of the spatial extended system taken in the 
same moments of time t\ and t<i is also small. The distance 
5(Ui,U2) has been defined as 



1/2 



S(Ui,U 2 ) 



|Ui(M)-U 2 (a, 



dx 



(6) 



where 1 1 • 1 1 is Euclidian norm. 

According to the neighbour method it means that there 
is the one-to-one correspondence between U(x,t) and y(t), 
therefore we can use the constructed auxiliary low di- 
mensional system y(t) to find UPTSTs by means of SD- 
method. 
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Fig. 1: The spatio-temporal dynamics of the charge density 
p(x,t) of the electron beam of Pierce diode. The oscillations 
for the selected control parameter value a = 2.8587T are chaotic 
both in space and time. 



Having constructed the auxiliary flow system (5) we 
can use SD-method to detect UPOs in it and UPSTSs 
in the initial spatial extended chaotic system (l)-(2), re- 
spectively. In R 3 space a plane p(x = 0.25, t) = 1.0 has 
been selected as Poincare secant. Let us denote the vec- 
tors y(t n ) = (1, p(0.5, t n ), p(0.75, t n )) T corresponding to 
the n-th crossing the selected secant surface by the tra- 
jectory y(t) as y n . Then the description of the system 
dynamics can be made with the help of the discrete map 



y n+ i = G(y n ), 



(7) 



where G(-) is the evolution operator. Obviously, it is im- 
possible to find the analytical form for the operator G, but 
numerical integration of the initial system of partial dif- 
ferential equations (l)-(2) can give us a sequence of values 
{y} n , generated by the map (7). 

SD-method for picking out unstable periodic orbits in 
the map (7) supposes consideration of the following map 
[23]: 



Yn+i = y„ + AC [G(y„) - y„] . 



(8) 



where A = 0.1 is the method constant and C is a certain 
matrix of the set C^. Each of matrices Ck should have only 
one non- vanishing entry +1 or —1 in row and column, i.e., 
they are orthogonal. In two dimensions the complete set 
of matrices consists of eight ones. 

In works [22,23] it was shown that map (8) under the 
appropriate choice of the matrix C allows to stabilize effec- 
tively the unstable saddle periodical orbits of systems (7) 
and (5). The positions of the UPOs in phase space are 
the same for the original chaotic system (7) and the trans- 
formed dynamical system (8) but their stability properties 
have changed: unstable fixed points turned into stable 
ones. A trajectory of transformed system (8) starting in 
the domain of attraction of a stabilized fixed point con- 
verges to it. Therefore, the UPOs of a chaotic dynamical 
system (7) can be obtained by iterating the transformed 
systems (8) using a robust set of initial conditions. In our 
calculation the matrix 



C = 



1 
1 




Fig. 2: The dependence of the space charge density p n (x — 
0.75) taken in the moments of time when the trajectory y(t) 
in R 3 space crosses the Poincare secant upon the number of 
iteration of SD-method for the UPSTS of the length 1 (T = 
4.2). Pierce parameter has been selected as a = 2.8587T. 



is suitable to find UPOs in (7). Having obtained UPOs for 
auxiliary systems (5) and (7) we can also obtain UPSTSs 
corresponding to them in the original spatial extended sys- 
tem (l)-(2). 

The transformed system (8) allows to find only the un- 
stable periodic orbits of length 1. To consider UPOs of 
length p the map 



Yn+i = y n + AC [G p (y n ) - y n ] 



(10) 



(9) 



should be considered instead of (8) where G^(-) is p- 
times iterated map (7). As far as the spatial extended 
system and the auxiliary flow system are considered, only 
the p-th crossing of the Poincare secant by the trajectory 
y(t) should be taken into account. 

So, by numerical iteration of the map (10) with differ- 
ent values of p one can find the set of the unstable peri- 
odic spatio-temporal states of the extended system (1)- 
(2). However, there is a problem concerning with search- 
ing the state XJ(x,t n +i) at the moment t n +i based on the 
known vector y n +i- Indeed, we know only the coordinates 
of the state y(t n+ i) in the Poincare secant but we don't 
know the corresponding distribution of space charge den- 
sity p(x, tn+i)j velocity v(x, t n +i) of the electron beam and 
the potential <p(x,£ n +i), and, correspondingly, we do not 
know the state U(x,£ n +i) of the extended system (1)- 
(2). However, as we have determined above with the 
help of the nearest neighbours method the state y(t n+ i) 
in the Poincare secant uniquely defines the correspond- 
ing spatial state U(x,t n+ i) belonging to the attracting 
finite-dimensional subspace W s of the infite-dimensional 
phase space W°° . To obtain this spatial state XJ(x,t n +i) 
mentioned above we have used the following procedure. 
The system of partial differential equations (l)-(2) de- 
scribing the fluid model of Pierce diode is integrated (and 
vector y(t) is calculated) untill some vector y(t s ) is close 
to the required one y n +i with some demanded precision: 
||y n +i — y( s )\\ < ^ where S is taken as 5 = 10 -3 . When 
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Fig. 3: The distribution of space charge density p(x,t) corre- 
sponding to the unstable spatio-temporal states with the fol- 
lowing lengthes p periods T: (a) p = 1, T = 4.2; (b) p = 2, 
T = 8.3; (c) p = 3, T = 16.9; (d) p = 4, T = 18.9. Pierce 
parameter has been selected as a — 2.8587T 



this condition is satisfied, the space state XJ(x,t s ) corre- 
sponding to the found vector y(s) are considered as the 
required one U(x,t n+ i) and then the next iteration ac- 
cording to (10) should be done. 

The spatio-temporal chaotic dynamics of the charge 
density p(x,t) of the electron beam of Pierce diode is 
shown in Fig. 1 for the Pierce parameter value a = 2.8587T. 
Applying the modified SD-method to the spatial extended 
system allows to find the demanded periodical time-space 
states. 

The convergence of the iteration procedure (10) is illus- 
trated by Fig. 2, which shows the dependence of the space 
charge density p n (x = 0.75) in the moments of time when 
the trajectory y(t) in R 3 space crosses the Poincare secant 
upon the number of iteration n of the SD-method when 
the unstable periodic spatio-temporal state of the length 
p = 1 is studied. One can see clearly that the iteration pro- 
cess of SD-method converges to the value corresponding to 
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Fig. 4: The spatio-temporal dynamics of the Ginzburg-Landau 
equation. The evolution of \u(x,t)\ is shown for the system 
length (a) L = 12.63 and (b) L = 13.25 



the unstable time-periodical spatio-temporal state of the 
system. Fig. 3 shows the distribution of space charge den- 
sity p(x, t) corresponding to the unstable spatio-temporal 
states with different periods T detected by means of SD- 
method. 

To verify both the correctness of the chosen value of the 
dimension m = 3 of the auxiliary system and the obtained 
results we have repeated the SD-method procedure for the 
value of the auxiliary system dimension being equal to 
m = 4. In this case the fourth-dimensional vector of the 
auxiliary is 



y(f) = (p(0.2,f),p(0.4,f),p(0.6,f), p(0.8, t)) 1 



(11) 



For m = 4 all UPSTSs founded coincide with the ones 
obtained above for the dimension of the auxiliary system 
m = 3, although the time of calculations increases in this 
case sufficiently. 

To show the universality of the proposed approach we 
also report the results of detecting the unstable periodic 
spatio-temporal states for the one-dimensional complex 
Ginzburg-Landau equation (CGLE) [32]. The CGLE is 
a fundamental model for the pattern formation and tur- 
bulence description. This equations is used frequently 
to describe many different nonlinear phenomena in laser 
physics [33] , chemical turbulence [34] , fluid dynamics [35] , 
bluff body wakes [36], coupled spatial extended systems 
[37,38]. 

We have considered one-dimensional CGLE 



du 



= u 



(l-ia)\u\ 2 u+ (l + i(3) 



(fiu 
dx 2 



(12) 



with periodical boundary conditions u(L,t) = u(0,t). All 
calculations were performed for a fixed system parameters 
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Fig. 5: The dependence of the value Reu n (x = L) upon the 
number of iteration of SD-method for UPSTS of the length 
p = 3 (T = 20.2). The values of Keu n (x = L) are taken 
in the moments of time when the trajectory y(t) in R 4 space 
crosses the Poincare secant. System length has been selected 
as L = 13.25. 




Fig. 6: The evolution of the module \u(x,t)\ corresponding to 
the unstable spatio-temporal states with the following lengths 
p and periods T: (a) p = 1, T = 12.1; (b) p = 2, T = 18.1. 
System length has been selected as L — 12.63. The dimension 
of the auxiliary vector y(t) has been chosen as m — 3. 



a = ft = 4 and random initial conditions. The numerical 
code was based on a semi-implicit scheme in time with 
finite differences in space. In all simulations we used a 
time step At = 0.0002 for the integration and a space 
discretization Ax = 0.04. 

The system length L has been chosen as the control 
parameter. In our study we examined two values of the 
control parameter: L\ = 12.63 and L2 = 13.25. For both 
these values of the control parameter L CGLE demon- 
strates the spatiotemporal chaotic regime. The corre- 
sponding spatio-temporal chaotic dynamics of CGLE are 
shown in Fig. 4 for the system lengths L = 12.63 and 
L = 13.25. One can see easily that the second case is 
characterized by more complex irregular spatio-temporal 
chaotic dynamics. Indeed, in the first case (L = 12.63) 
the chaotic dynamics is characterized by only one positive 
Lyapunov exponent Ai = 0.04, while the second chaotic 
regime (L = 12.63) is characterized by two positive Lya- 
punov exponents Ai = 0.10 and A2 = 0.07. 

Applying the modified SD-method to the spatial ex- 
tended CGLE we can find the demanded unstable period- 
ical spatio-temporal states as well as for the fluid model 
of Pierce diode. We have constructed the vector (5) of the 
auxiliary low dimensional system as 

y(t) = (u(x 1 ,t),...,u(x m ,t)) T , (13) 

where m is the dimension of the auxiliary system vector, 
Xi = iL/m, i = 1, m. 

In contrast to the fluid model of Pierce diode (l)-(2) 
the dimension m of the auxiliary vector y(t) is unknown 
for CGLE (12). Therefore, we have to try to find UPSTSs 
by means of the SD-method (10) for the different values of 
the auxiliary system dimension m starting from the mini- 
mal dimension value m = 3. If the required UPSTS is not 
found for the selected value of the auxiliary vector dimen- 
sion m*, the SD-method procedure should be repeated for 



the greater dimension value m = m* + 1. 

For the system length L = 12.63 the dimension of the 
auxiliary system m = 3 is found to be adequate for the 
correct UPSTSs detection. As it was mentioned above the 
system behavior is characterized by one positive Lyapunov 
exponent. For the more complicated case L = 13.25 (when 
the behavior of CGLE is characterized by two positive 
Lyapunov exponents) the dimension of the auxiliary vector 
should be taken as m = 4 for UPSTSs to be detected 
successfully. The matrixes 

C=("J J) and C=^0 1 j (14) 

are found to be suitable to find UPSTSs for the system 
lengths L = 12.63 and L = 13.25, respectively. 

The convergence of the iteration procedure (10) is illus- 
trated in Fig. 5. One can see clearly that the iteration pro- 
cess of SD-method converges to the value corresponding to 
the unstable periodic spatio-temporal state of the system. 
Fig. 6 shows the evolution of the profiles \u(x,t)\ corre- 
sponding to the unstable periodic spatio-temporal states 
with the different periods T detected by means of SD- 
method for the system length L = 12.63, when the di- 
mension of the auxiliary vector (13) has been chosen as 
n = 3. The analogous evolution of the profiles \u(x, t)\ cor- 
responding to the unstable periodic spatio-temporal states 
with the different lengths p and periods T is shown in 
Fig. 7 for L = 13.25 and m = 4. 

In conclusion, we have proposed the method of the de- 
tection of the unstable periodic spatio-temporal states of 
spatial extended chaotic systems being the extension of the 
well known SD-method. The effectiveness of this method 
is illustrated by the consideration of the fluid model of 
Pierce diode and the complex Ginzburg-Landau equation. 
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Fig. 7: The evolution of the proviles \u(x,t)\ corresponding to 
the unstable periodic spatio-temporal states with the following 
lengths p and periods T: (a) p = 1, T = 4.95; (b) p = 2, 
T = 16.83; (c) p = 3, T = 20.20. The system length has been 
selected as L = 13.25. The dimension of the auxiliary vector 
y(t) has been chosen as m — 4 
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